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We present a search for the standard model (SM) Higgs boson produced in association with 
a Z boson in 9.7 fb _1 of pp collisions collected with the DO detector at the Fermilab Tevatron 
Collider at yfs = 1.96 TeV. Selected events contain one reconstructed Z — ¥ e + e~ or Z — > fi + n~ 
candidate and at least two jets, including at least one jet identified as likely to contain a b quark. 
To validate the search procedure, we also measure the cross section for ZZ production in the same 
final state. It is found to be consistent with its SM prediction. We set upper limits on the ZH 
production cross section times branching ratio for H — > bb at the 95% C.L. for Higgs boson masses 
90 < M H < 150 GeV. The observed (expected) limit for M H = 125 GeV is 7.1 (5.1) times the SM 
cross section. 

PACS numbers: 13.85.Ni, 13.85.Qk, 13.85.Rm, 14.80.Bn 



In the standard model (SM), the spontaneous break- 
ing of the electroweak gauge symmetry generates masses 
for the W and Z bosons and produces a residual mas- 
sive particle, the Higgs boson Precision electroweak 
data, including the latest W boson mass measurements 
from the CDF Q and DO Q Collaborations, and the 
latest Tevatron combination for the top quark mass [4| 
constrain the mass of the SM Higgs boson to Mjj < 
152 GeV 0] at the 95% confidence level (CIA Direct 
searches at the CERN e+e" Collider (LEP) [6], by the 
CDF and DO Collaborations at the Fermilab Tevatron pp 
Collider Q , and by the ATLAS and CMS Collaborations 
at the CERN Large Hadron Collider (LHC) d©] further 
restrict the allowed range to 116.6 < Mjj < 119.4 GeV 
and 122.1 < M H < 127.0 GeV. The ATLAS and CMS 
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results indicate excesses above background expectations 
at Mjj ~ 125 GeV. With additional data and analysis 
improvements, the LHC experiments confirm their initial 
indications and observe a particle with properties consis- 
tent with those predicted for the SM Higgs boson [lfj. 

For Mh 5, 135 GeV, the primary decay is to the bb 
final state [1JJ. At the Tevatron, the best sensitivity to 
a SM Higgs boson in this mass range is obtained from 
the analysis of its production in association with a W 
or Z boson and its subsequent decay into bb. Evidence 
for a signal in this decay mode would complement the 
LHC findings and provide further indication that the new 
particle is the SM Higgs boson. 

We present a search for ZH — > £ + £~bb events, where I 
is cither a muon or an electron. The data for this analysis 
were collected at the Tevatron at y/s = 1.96 TeV with 
the DO detector from April 2002 to September 2011 and 
correspond to an integrated luminosity of 9.7 fb _1 after 
data quality requirements are imposed, which represents 
the full Run II data set. To validate the search proce- 
dure, we also present a measurement of the ZZ produc- 
tion cross section in the same final states and topologies 
used for the search. The results presented here supersede 
our previous search in the ZH — ¥ £ + £~bb channel [l2j]. 
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Beyond the inclusion of additional data, the most signif- 
icant updates to this analysis are the use of an improved 
6-jet identification algorithm, revisions to the kinematic 
fit, and a new multivariate analysis strategy. A search 
for ZH — > £ + £~bb has also been performed by the CDF 
Collaboration [l3j . 

The DO detector 0, consists of a central tracking 
system within a 2 T superconducting solenoidal magnet 
and surrounded by a preshower detector, three liquid- 
argon sampling calorimeters, and a muon spectrometer 
with a 1.8 T iron toroidal magnet. In the intercryostat 
regions (ICRs) between the central and end calorime- 
ter cryostats, plastic scintillator detectors enhance the 
calorimeter coverage. The analyzed events were ac- 
quired predominantly with triggers that select electron 
and muon candidates online. However, events satisfying 
any trigger requirement are considered in this analysis. 

The event selection requires a pp interaction vertex 
that has at least three associated tracks. Selected events 
must contain a Z — > £ + £~ candidate. The analysis is 
conducted in four separate channels. The dimuon (h/j,) 
and dielectron (ee) channels include events with at least 
two fully reconstructed muons or electrons. In addition, 
muon-plus-track (/x/Ztrk) and electron-plus-ICR electron 
(eeicR.) channels are designed to recover events in which 
one of the leptons points to a poorly instrumented region 
of the detector. 

The event selection requires at least two muons 
identified in the muon system, both matched to central 
tracks with transverse momenta pr > 10 GeV. At least 
one muon must have \r)\ < 1.5, where rj is the pseudo- 
rapidity, and pt > 15 GeV. At least one of the muons 
must be separated from any jet with px > 20 GeV and 
\r)\ < 2.5 by ATI = ^ Arj 2 + Ac/) 2 > 0.5, from other 
tracks, and from energy deposited in the calorimeter. We 
also apply isolation requirements based on the ratios of 
the calorimeter energy and the sum of pt of tracks near 
the lepton to the lepton px in this analysis. 

The /x/itrk event selection requires exactly one muon 
with 1 77 1 < 1.5 and pr > 15 GeV that is isolated both in 
the tracker and in the calorimeter. In addition, a second 
isolated track reconstructed in the tracker with \r]\ < 2 
and pr > 20 GeV must be present. Its distance A1Z from 
the muon and from any jet of pt > 15 GeV and \r)\ < 2.5 
must be greater than 0.1 and 0.5, respectively. For the 
/i/j. and /i/Jtrk channels, the two muon-associatcd tracks 
must have opposite charge. 

The ee event selection requires at least two electrons 
with transverse energy Et > 15 GeV that pass selec- 
tion requirements based on the energy deposition and 
shower shape in the calorimeter and the preshower de- 
tector. Both electrons are required to be isolated in the 
tracker and the calorimeter. At least one electron must 
be identified in the region |^| < 1.1. The electrons in 
< 1.1 must match central tracks or a set of hits in the 
tracker consistent with that of an electron trajectory. 

The eeicR event selection requires exactly one electron 
in the calorimeter with Et > 15 GeV and a track point- 



ing toward one of the ICRs, 1.1 < < 1.5. The track 
must be isolated, be matched to a calorimeter energy de- 
posit with Et > 10 GeV and have px > 15 GeV. For 
the ee and eeicR selections, electrons must be separated 
from all jets by ATI > 0.5. 

Jets are reconstructed in the calorimeter by using the 
iterative midpoint cone algorithm [1 6j j with a cone of ra- 
dius 0.5 in rapidity and azimuthal angle. The jet identi- 
fication efficiency is « 95% at pt = 20 GeV and reaches 
99% at pt = 50 GeV. Jets are denoted as "taggable" 
if the associated tracks meet criteria that algorithms to 
identify jets as likely to contain 6-quarks operate effi- 
ciently. The taggability efficiency is at least 90% for most 
of the jets in this analysis. We use "inclusive" to denote 
the event sample selected by requiring the presence of 
two leptons and use "pretag" for the event sample that 
meets the additional requirements of having at least two 
taggable jets with pt > 20 GeV and |?7| < 2.5 and a 
dilepton invariant mass 70 < mu < 110 GeV (l7j . 

Jets are identified as likely to contain b quarks (6- 
tagged) if they pass "loose" or "tight" requirements on 
the output of a multivariate discriminant trained to sep- 
arate b jets from light jets. This discriminant is an im- 
proved version of the neural network ^-tagging discrimi- 
nant described in Ref. [lsj . For taggable jets in \rj\ < 1.1 
and with pt ~ 50 GeV, the 6-tagging efficiency for b 
jets and the misidentification probability of light (uds or 
gluon) jets are, respectively, 72% and 6.7% for loose b 
tags, and 47% and 0.4% for tight b tags. Events with 
at least one tight and one loose b tag are classified as 
double-tagged (DT). Events not in the DT sample that 
contain a single tight b tag are classified as single-tagged 
(ST). 

The dominant background process is the production 
of a Z boson in association with jets, with the Z decay- 
ing to dileptons (Z+jets). The light-flavor component 
(Z+LF) includes jets from only light quarks or gluons. 
The heavy-flavor component (Z+HF) includes Z + bb, 
which has the same final state as the signal, and Z + cc 
production. The remaining backgrounds are from ti pro- 
duction; WW, WZ, and ZZ (diboson) production; and 
multijet (MJ) events with nonprompt muons or with jets 
misidentified as electrons. 

We simulate ZH and diboson production with 
PYTHIA [l9[ . In the ZH samples, we consider the con- 
tributions to the signal from the £ + £~bb, £ + £~cc, and 
£+£-t+t- final states. The £ + £~bb accounts for 99% 
(97%) of thesignal yield in the DT (ST) sample. The 
Z+jets and tt processes are simulated with ALPGEN [2(| , 
followed by PYTHIA for parton showering and hadroniza- 
tion [2l|. All simulated samples are generated by us- 
ing the CTEQ6L1 22] leading-order parton distribution 
functions. We process all samples by using a detector 
simulation program based on GEANT3 [23| and the same 
offline reconstruction algorithms used for data. We over- 
lay events from randomly chosen beam crossings with 
the same instantaneous luminosity distribution as data 
on the generated events to model the effects of multiple 
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pp interactions and detector noise. 

We take the cross sections and branching ratios for 
signal from Refs. [ll|, H3]- For the diboson processes, 
we use next-to-leading order (NLO) cross sections from 
the Monte Carlo program MCFM |25|. We scale the it 
cross section to approximate next-to-NLO (2|| and the 
inclusive Z boson cross section to next-to-NLO [27j and 
apply additional NLO heavy-flavor corrections to the Z+ 
bb and Z + cc samples, calculated from mcfm to be 1.52 
and 1.67, respectively. 

To improve the modeling of the pt distribution of the 
Z boson, we reweight simulated Z+jets events to be con- 
sistent with the measured pr spectrum of Z bosons in 
the data [Hj]. We correct the energies of simulated jets 
to reproduce the resolution and energy scale observed in 
the data [29j]. We apply the trigger efficiencies, measured 
in the data, as event weights to the simulated p.p, 1 /x/Xtrk 
and eeicR events. In the ee channel, we have verified that 
the trigger efficiency is consistent with 100% for our se- 
lection. We apply scale factors to account for differences 
in reconstruction efficiency between the data and simula- 
tion. Motivated by a comparison with data [30( and the 
SHERPA generator [3l|, we reweight the Z+jets events to 
improve the ALPGEN modeling of the distributions of the 
r\ of the two jets. 

We estimate the M J backgrounds from control samples 
in data obtained by inverting some of the lepton selec- 
tion requirements, e.g., the lepton isolation requirements 
in the /i/i channel and the shower shape requirements 
in the ee channel. We adjust the normalizations of the 
MJ background and all simulated samples by scale fac- 
tors determined from a simultaneous fit to the m« dis- 
tributions in the 0-jet, 1-jet, and > 2-jet samples of each 
lepton selection. The inclusive sample constrains the lep- 
ton trigger and identification efficiencies, while the pre- 
tag sample, which includes jet requirements, is used to 
correct the Z+jets cross section. The total event yields 
after applying all corrections and normalization factors 
are shown in Table |TJ The observed event yields are con- 
sistent with the expected background. 

To exploit the fully constrained kinematics of the 
ZH — > £ + £~bb process, we adjust the energies of the can- 
didate leptons and jets within their experimental resolu- 
tions by using a likelihood fit that constrains ma to the 
mass and width of the Z boson and constrains the pr of 
the £ + £~bb system to zero with an expected width deter- 
mined from ZH Monte Carlo events. This kinematic fit 
improves the dijet mass resolution by 10%— 15%, depend- 
ing on Mh ■ The dijet mass resolution for Mh — 125 GeV 
is 15 GeV with the kinematic fit [171 ] . 

We use a two step multivariate analysis strategy based 
on random forest (RF, an ensemble classifier that con- 
sists of many decision trees) discriminants [32^ . as im- 
plemented in the TMVA software package 33 1 , to improve 
the separation of signal from background 17J • We choose 
well modeled kinematic variables that are sensitive to the 
ZH signal as inputs for the analysis. These include the 
Pt of the two &-jet candidates and the dijet mass, before 



and after the jet energies are adjusted by the kinematic 
fit. In the first step, we train a dedicated RF (it RF) that 
takes it as the only background and ZH as the signal. 
This approach takes advantage of the characteristic sig- 
nature of the it background, for instance, the presence 
of large missing transverse energy. In the second step, 
we use the it RF to define two independent regions: a it 
enriched region (it RF < 0.5) and a it depleted region (it 
RF > 0.5). The it depleted region contains 94% (93%) of 
the DT (ST) signal contribution and 55% (82%) of DT 
(ST) background events. In each region, we train a global 
RF to separate the ZH signal from all backgrounds. In 
both steps we consider ST and DT events separately and 
train the discriminants for each assumed value of Mh in 
5 GeV steps from 90 to 150 GeV. 

We assess systematic uncertainties resulting from the 
background normalization for the MJ contribution, typ- 
ically 10%. The normalization of the Z+jets sample to 
the pretag data constrains that sample to the statisti- 
cal uncertainty, <1%, of the pretag data. Because this 
sample is dominated by the Z+LF background, the nor- 
malization of the it, diboson, and ZH samples acquires 
a sensitivity to the inclusive Z cross section, for which 
we assess a 6% uncertainty [27|. We assign this uncer- 
tainty to these samples as a common uncertainty. For 
Z+HF, a cross section uncertainty of 20% is determined 
from Ref. [25| . For other backgrounds, the uncertain- 
ties are 6%-10% [25|, [26l. For the signal, the cross sec- 
tion uncertainty is 6% [24|. Sources of systematic un- 
certainty affecting the shapes of the final discriminant 
distributions are the jet energy scale, l%-3%; jet energy 
resolution, 2%-4%; jet identification efficiency, w 4%; 
and 6-tagging efficiency, 4%-6%. Other sources include 
trigger efficiency, 4%-6%; parton distribution function 
uncertainties [341 ] . <1%; data-determined corrections to 
the model for Z+jets, 3%-4%; modeling of the underly- 
ing event, <1%; and from varying the factorization and 
renormalization scales for the Z+jets simulation, <1%. 

The global RF distributions from the four samples (ST 
and DT in the it depleted and it enriched regions) in each 
channel along with the corresponding systematic uncer- 
tainties are used for the statistical analysis of the data. 
We set 95% C.L. upper limits on the ZH cross section 
times branching ratio for H — > bb with a modified frc- 
quentist (CL S ) method that uses the log likelihood ra- 
tio of the signal+background (S+B) hypothesis to the 
background-only (B) hypothesis 35]. To minimize the 
effect of systematic uncertainties, we maximize the likeli- 
hoods of the B and S+B hypotheses by independent fits 
that allow the sources of systematic uncertainty to vary 
within their Gaussian priors [36| . 

To validate the search procedure, we search for ZZ 
production in the £ + £~bb and £ + £~cc final states. We 
use the same event selection, corrections to our signal 
and background models, and RF training procedure as 
for the ZH search 17]. Our search also includes WZ 
production in the cs£ + £~ final state. We collectively re- 
fer to these as V Z production. Using the same modified 
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TABLE I: Expected and observed event yields for all lepton channels combined after requiring two leptons (inclusive), after 
also requiring at least two jets (pretag), and after requiring exactly one (ST) or at least two (DT) 6-tags. The ZH signal yields 
are for Mh = 125 GeV. The uncertainties quoted on the total background for ST and DT and signal include the statistical 
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MJ 
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17.3 ± 
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9.2 ± 
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824 ± 102 


54 


60 


600 


33 
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366 ± 39 
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FIG. 1: (color online). Distributions of the global RF discrim- 
inant in the tt depleted region, assuming Mh = 125 GeV, 
after the fit to the background-only model for data (points 
with statistical error bars) and background (histograms) for 
(a) single-tagged events and (b) double-tagged events, (c) 
Background-subtracted distribution for (b). The signal dis- 
tribution is shown with the SM cross section scaled by a factor 
of five. The blue lines indicate the uncertainty from the fit. 



frequentist method as for the ZH search and fitting the 
RF distributions to the S+B hypothesis, we measure a 
VZ cross section of 0.8 ± 0.4 (stat) ± 0.4 (syst) times that 
of the SM prediction with a significance of 1.5 standard 
deviation (s.d.) and an expected significance of 1.9 s.d. 
This result is consistent with the recent DO ZZ + WZ 
cross section measurement obtained in fully leptonic de- 



cay channels [37J . 

The output of the RF trained to separate signal events 
with Mh = 125 GeV from background is shown in Fig.Q] 
for ST and DT events separately in the tt depleted re- 
gion, after the background-only fit. Also shown is the 
background-subtracted RF distribution for DT events in 
the data. The upper limit on the cross section times 
branching ratio for H — > bb, expressed as a ratio to the 
SM prediction, is presented as a function of Mg in Ta- 
ble [CD and Fig. [1 At M H = 125 GeV, the observed 
(expected) limit on this ratio is 7.1 (5.1). The expected 
limits are « 20% lower than those anticipated from the 
increase in the data because of the analysis improvements 
described above. 



■ Observed Limit 
Expected Limit 

1 Expected ± 1 s.d 
Expected ± 2 s.d 



9.7 fb" 




l-Q 
-Q 

CD 

X 10 

S 

,1 



130 140 150 

M H (GeV) 

FIG. 2: (color online). Expected and observed 95% C.L. cross 
section upper limits on the ZH cross section times branching 
ratio for H — > bb, expressed as a ratio to the SM prediction. 



In summary, we have searched for SM Higgs boson 
production in association with a Z boson in the final state 
of two charged leptons (electrons or muons) and two b- 
quark jets by using a 9.7 fb -1 data set of pp collisions at 
y/s = 1.96 TeV. We also measure the cross section for 
VZ production in the same final state with the result of 
0.8 ± 0.4 (stat) ± 0.4 (syst) times its SM prediction. We 
set an upper limit on the ZH production cross section 
times branching ratio for H — > bb as a function of Mh- 
The observed (expected) limit for Mh = 125 GeV is 7.1 
(5.1) times the SM cross section. 
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TABLE II: The expected and observed 95% C.L. upper limits on the ZH production cross section times branching ratio for 
ZH — > l + £~bb, expressed as a ratio to the SM prediction. 



M H (GeV) 
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Expected 


2.6 


2.7 


2.8 


3.0 


3.4 
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12 


18 


29 
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1.8 


2.3 


2.2 
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6.2 
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12 


16 


19 


31 
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Appendix 
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FIG. 3: The dilepton mass spectra in the (a) n\x, (b) /x/Jtrk, (c) ee and (d) eeicR channels. Distributions are shown in the 
pretag control sample, in which all selection requirements except 6-tagging are applied. Signal distributions, for Mh = 125 
GeV, are scaled by a factor of 500. 

The dimuon and dielectron mass spectra, after requiring two leptons and at least two jets are shown in Fig. [3] 
Distributions of the dijet invariant mass spectra before and after adjustment by the kinematic fit, are shown in Fig. 2) 
A complete list of RF input variables is shown in Table IIIII Comparisons of the data and MC distributions of the 
ti RF output summed over all lepton channels are shown for Mh — 125 GeV in Figure [5] Post-kinematic fit dijet 
mass distributions for ST and DT in the tt depleted region are shown in Fig. |U Fig. [7] displays the global RF 
distributions in the ti enriched region, after the fit to the background-only hypothesis. Fig. [5] shows the observed 
LLR as a function of Higgs boson mass. Also shown are the expected (median) LLRs for the background-only and 
signal+background hypotheses, together with the one and two standard deviation bands about the background-only 
expectation. Fig. ® shows the post-fit RF distributions in the tt depleted region for the VZ search. Fig. [TU] displays 
the post-fit distribution of the dijet invariant mass from the kinematic fit. 
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FIG. 4: Dijet invariant mass distributions before the kinematic fit in (a) ST events and (b) DT events; and after the kinematic 
fit in (c) ST events and (d) DT events, combined for all lepton channels. Signal distributions (Mb = 125 GeV) are shown with 
the SM cross section multiplied by 20. 



TABLE III: Variables used for the tt and global RF training. The jets that form the Higgs boson candidate are referred to as 
b\ and 62. 



variables 


definition 


tt RF 


global RF 


m bb (.m bb lxt ) 


invariant mass of the dijet system before (after) the kinematic fit 


V 


V 


Pt (Pt ) 


transverse momentum of the first jet before (after) kinematic fit 




V 


Pt (Pt J 


transverse momentum of the second jet before (after) kinematic fit 


V 


V 


n bb 

Pt 


transverse momentum of the dijet system before the kinematic fit 


V 


V 




A<f) between the two jets in the dijet system 




V 


Ar?(6i,6 2 ) 


Arj between the two jets in the dijet system 








invariant mass of all jets in the event (the multijet mass) 


V 


V 




transverse momentum of all jets in the event 




V 




scalar sum of the transverse momenta of all jets in the event 


V 




Pt/(\Pt\ + \Pt\) 


ratio of dijet system pt over the scalar sum of the pt of the two jets 


V 




m u 


invariant mass of the dilepton system 


V 




Pt 


transverse momentum of the dilepton system 


V 


V 




A(j) between the two leptons 


V 


V 


colinearity(£i,f2 


cosine of the angle between the two leptons (colinearity) 


V 


V 


a^u, 66) 


A<j> between the dilepton and dijet systems 




V 


cos#* 


cosine of the angle between the incoming proton and the Z in the zero momentum frame a 




V 


m(Ubb) 


Invariant mass of dilepton plus dijet system 






H T (Ubb) 


Scalar sum of the transverse momenta of the leptons and jets 




V 


$T 


missing transverse energy of the event 


V 




Pt 


the lf,T significance b 


V 


V 


— In Lfu 


negative log likelihood from the kinematic fit 


V 




tt RF 


tt RF output 




V 



a S. Parke and S. Veseli, Phys. Rev. D 60, 093003 (1999). 
b A. Schwartzman, FERMILAB-THESIS-2004-21. 
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tt RF Output tt RF Output 

FIG. 5: The tt RF output (Mh = 125 GeV) for all lepton channels combined (a) ST and (b) DT events. Signal distributions 
are shown with the SM cross section multiplied by 20. The vertical arrows indicate the tt RF selection requirement used to 
define the tt enriched and depleted samples. 




FIG. 6: Post-kinematic fit dijet mass distributions in the tt depleted region for all lepton channels combined assuming Mh = 125 
GeV for (a) ST events and (b) DT events. Signal distributions are shown with the SM cross section multiplied by 20. 
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FIG. 7: Post-fit RF output distributions in the tt enriched region, assuming Mb = 125 GeV, after the fit to the background-only 
model for (a) ST events and (b) DT events. Background-subtracted distributions for (a) and (b) are shown in (c) and (d), 
respectively. Signal distributions are shown with the SM cross section scaled to 50 x SM prediction in (c) and (d). The blue 
lines are the total posterior systematic uncertainty following a fit of the background to the data. 
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FIG. 8: Observed LLR as a function of Higgs boson mass. Also shown are the expected LLRs for the background-only 
(B) and signal+background (S+B) hypotheses, together with the one and two standard deviation (s.d.) bands about the 
background-only expectation. 
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FIG. 9: Post-fit VZ RF output distributions in the it depleted region after the fit to the S+B model for (a) ST events and (b) 
DT events. Background-subtracted distributions for (a) and (b) are shown in (c) and (d), respectively. Signal distributions are 
scaled to the measured VZ cross section. The blue lines indicate the uncertainty from the fit. 
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FIG. 10: Post-fit distributions of the dijet invariant mass (from the kinematic fit) in the tt depleted region after the fit to the 
S+B model for (a) ST events and (b) DT events. Background-subtracted distributions for (a) and (b) are shown in (c) and (d), 
respectively. Signal distributions are scaled to the measured VZ cross section. The blue lines indicate the uncertainty from 
the fit. 



